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We  have  studied  the  parametric  instabilities  of  electron 
cyclotron  waves  in the  ^BT^ring  and^JJ)  the  large  tokamaks 
eg.  PLT.  In  the  EBT,  the  electron  cyclotron  pump  of  finite 
wavenumber  k^  decays  into  two  Bernstein  modes  at  the  second 
harmonic  cyclotron  layer  and  can  account  for  the  heating  of 
the  ring  in  the  initial  phase.  The  coupling  coefficient 
for  this  decay  vanishes  for  a  dipole  pump,  whereas  the 
convective  threshold  with  finite  k^  is  ^200  W/cjh^.  For  large 
tokamaks,  the  convective  threshold  for  various  decay  channels 


2 


turns  out  to  be  >200  KW/cw  at  3ram  wavelength.. 
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I .  INTRODUCTION 


The  recent  developments  of  powerful  microwave  sources  in  the  millimeter 
range  have  aroused  substantial  efforts  on  the  electron  cyclotron  heating  of 
large  devices,  eg.,  Elmo  bumpy  torus,  tokamak  and  mirrors.  In  the  Elmo 
bumpy  torus^  (EBT)  the  microwaves  are  substantially  absorbed  at  the  second 
harmonic  cyclotron  layer  (ojo  =  2^)  and  a  hot  electron  ring  with  several 
keV  temperature  is  formed.  No  satisfactory  explanation  of  how  the  ring 
is  formed  is  available  so  far.  Here  we  examine  a  possible  parametric  in¬ 
stability  at  the  second  harmonic  which  may  play  a  role  in  the  acceleration 

of  the  high  energy  electrons. 
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In  the  case  of  tokamak  ’  ,  the  experiments  on  electron  cyclotron 
heating  of  ISXB,  TM-3,  Tuman-2  and  other  devices  have  proved  promising 
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with  almost  a  hundred  percent  power  absorption  efficiency.  Ott  et  al. 

and  Eldridge  et  al.^  have  predicted  the  absorption  efficienty  of  100%  for 

the  ordinary  and  extraordinary  modes  in  PLT  on  the  basis  of  linear  theory 

of  cyclotron  absorption.  The  presently  employed  power  densities  for  tokamak 
2 

heating  (<  1  KW/cm  )  at  electron  cyclotron  frequency  are  far  too  low  to 

excite  any  parametric  processes.  Nevertheless  the  use  of  higher  power 

dnesities  in  future  experiments  might  initiate  such  processes.  An  estimate 

of  these  effects  is  necessary  for  the  prediction  of  the  heating  rate. 
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Ott  et  al .  have  recently  studied  the  parametric  instabilities  of  the 
ordinary  pump  wave  in  PLT  on  the  low  magnetic  field  side  of  the  torus; 
the  parametric  instabilities  of  the  extraordinary  mode  have  been  omitted 
on  the  ground  that  on  the  inner  side  of  the  torus,  from  where  this  mode 
is  launched,  the  electrostatic  high  frequency  modes  possess  frequencies 
greater  than  that  of  the  pump.  The  two  channels  of  decay  viz,  resonant 
decay  into  upper  hybrid  and  ion  cyclotron  waves  and  nonresonant  decay  into 


electron  Bernstein  mode  have  been  studied.  The  threshold  powers  for  these 
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channels  are  greater  than  200  KW/cm  ,  hence  they  are  unimportant.  Elder 
and  Perkins^  have  investigated  the  parametric  decay  of  the  extraordinary 
mode,  launched  from  the  outside  of  the  torus,  into  ion  acoustic  and  Langmuir 
waves.  The  threshold  power  for  this  decay  is  very  large:  powers  of  the 
order  of  350  MW  are  required  for  its  onset  in  PDX  device. 

In  this  paper,  we  have  investigated  the  parametric  instabilities  of 
electron  cyclotron  waves  I)  at  the  second  harmonic  cyclotron  layer  in  EBT 
and  II)  at  the  cyclotron  harmonic  in  PLT.  In  the  EBT,  the  pump  wave  could 
decay  into  two  Bernstein  waves  which  possess  slow  group  velocities  and  the 
decay  channel  should  have  lowest  convective  threshold.  In  the  case  of  PLT, 
we  shall  match  the  frequency  of  the  pump  to  the  cyclotron  frequency  at  the 
center  of  the  torus  (x=0)  and  would  be  interested  in  the  parametric  in¬ 


stabilities  in  the  vicinity  of  the  center.  The  extraordinary  mode  launched 


from  the  inside  of  the  torus  (x=-rQ)  could  decay  into  i)  two  lower  hybrid 
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waves,  ii)  lower  hybrid  and  ion  cyclotron  waves  (for  <  2u>c)  and 

2  2  2 

iii)  lower  hybrid  and  quasimodes  (for  a>c  <  <  2a>c)  on  the  innerside  of  the 

torus.  The  ordinary  mode  launched  from  the  outside  of  the  torus  (x=rQ)  could 

decay  into  i)  two  lower  hybrid  waves,  ii)  Bernstein  and  quasimodes,  iii)  upper 

2  2 

hybrid  and  ion  cyclotron  modes  and  quasimodes  (for  <  uic)  .  The  last 
channel  is,  however,  possible  in  the  low  density  region  hence  we  would  ignore 


it. 


In  Section  II  we  have  obtained  the  coupling  coefficient  for  the  decay 


of  a  finite  wavenumber  (kQ}*0)  electron  cyclotron  pump  around  u)Q  «  2u>c  into 
two  Bernstein  modes  in  EBT  by  solving  the  Vlasov  equation  for  electrons  in 
the  guiding  center  coordinates.  The  coupling  coefficient  vanishes  for  ko”0. 
Following  Liu^,  the  convective  threshold  for  the  instability  has  been  ob¬ 
tained.  In  Section  III,  we  have  studied  the  decay  instability  of  electron 


r 
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cyclotron  pump  (tu0  -  wc)  in  a  tokamak  by  employing  fluid  theory,  which  is 
applicable  for  long  wavelength  decay  modes.  A  discussion  of  results  is 
given  in  Section  IV. 
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II.  TWO  BERNSTEIN  MODE  DECAY  IN  EBT 


We  consider  the  propagation  of  an  extraordinary  pump  wave 


£  exp  -  i((o  t  -  k  x) 
o  o  o 


into  the  mirror  midplane  of  EBT,  transverse  to  the  direction  of  static  magnetic 

field  S  II  z;  E  =  -iu  u2E  /oi  (u)2  -  u2),  E  =  0,  B  =  k  cE  /w  ,  B„  = 
s  11  ox  c  p  oy  o  o  c  *  oz  *  oz  o  oy  o’  ox 
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o  c  1  2  2  2  2J  K  1 
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and  oip  and  oi c  are  the  electron  plasma  and  cyclotron  frequencies.  The  density 
gradient  is  along  x  axis  and  the  pump  frequency  is  close  to  2wc.  The  pump 
wave  decays  into  two  electrostatic  waves  E(w,k)  =  -V[<(>exp  -  i(ut  -  k  •  x)] 
and  E^w^k^)  =  -V^exp  -  i(o)^t  -  kj  •  x)  J  where  Wj  ■  w  -  u>o,  k^  *  k  -  kQ. 
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The  response  of  electrons  to  these  fields  is  governed  by  the  Vlasov  equation 


3f  ,  Pz  3f  •  3f  •  3f  _T  3f  .  •  3f  .  •  »f 

— —  +  —  — —  +  9  —  +  p  - —  -  eE  - —  +  x  -5 —  +  y  ~ - 

3 1  m  3z  30  3p  z  3p  g  3x  'g  3y 

z  K  8 


where  u  *  mvi/2u>c  is  the  magnetic  moment,  0  is  the  gyrophase  angle,  v^  is  the 

perpendicular  velocity,  -e  and  m  are  the  electronic  change  and  mass,  respectively, 

a)  *  eBs/mc  and  x  =  x  -  v,  x  oi  /to 2  are  the  guiding  center  coordinates.  Using 
c  g  ICC 

the  equation  of  motion,  one  can  write 
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where  the  superscript  T  refers  to  the  total  field.  Now  we  express  the  total 
electric  field  and  the  distribution  function  as 
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where  J  =  J  (k  P) ,  J°(k  P) ,  J  =  J  (k..  p),  p  =  v./u>  ,  prime  over  the 
ii  n  jl  n  o  n  x  J*  x  c 

Bessel  function  denotes  the  derivative  with  respect  to  its  argument,  6q  is 

~y  -> 

the  angle  between  Eq^  and  x  axis,  d  between  k^  and  x  axis  and  d^  between 
~y 

k.^  and  x  axis.  Now  onwards  we  shall  assume  that  k|(  =  0  for  the  Bernstein 
modes.  Using  Eqs.  (4)  -  (8)  in  Eq.  (3),  we  obtain  the  linear  response. 
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where  we  have  assumed  fQ  to  be  Maxwellian  at  temperature  T.  For  the  nonlinear 
response  we  get, 
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First,  let  us  mention  here  that  in  the  limit  of  kQ  =  0  (dipole  pump),  f^ 

NL 

is  greatly  simplified.  However,  the  nonlinear  density  perturbation  n 
obtained  from  this  after  carrying  out  the  integration,  vanishes,  identi¬ 
cally.  In  the  next  approximation  we  retain  terms  that  go  linearly  with  kQ 
and  possess  resonant  denominators,  then  we  obtain 


CO 
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where  nQ  is  the  unperturbed  electron  density,  b  =  ^ive^a)c’  ve  =  (^T/m) 
and  1^  is  the  modified  Bessel  function  of  argument  b.  Similarly,  one 
obtains, 

2.  V  .  'i<S  '  S°  +  V 

NL  o  6  oi  o  ^  6  d  _-b, 
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Using  Eqs.  (11)  and  (12)  in  the  Poisson's  equation,  we  get  the  following 
nonlinear  dispersion  relation 
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1  =  y/e  c1  ,  (13) 

where  v  Is  coupling  coefficient  and  c,  £l  are  the  dielectric  functions 
at  (w,k),  (u^.k.^); 

(14) 

(15) 

Using  the  linear  dispersion  relation  (e  =  0)  for  the  Bernstein  mode,  one 
may  take 
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and  hence 
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u  =  eEoi/mwc.  In  writing  Eq.  (16)  we  have  assumed  b  >  1.  Following 
Liu,  the  convective  amplification  of  the  decay  waves  from  Eq.  (13)  can  be 
obtained  as 
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Taking  k  v  /w  'v  3,  v  / u>  ^  10  2,  k2L  /k  'v  10^,  k  ■v  3  x  10  and  expressing 
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2 

p  in  KW/cm  ,  we  get 
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Thus,  the  threshold  power  for  two  Bernstein  mode  decay  is  =200  W/cm  for 
(toc/ujp)  =  2.  For  shorter  wavelength  modes  (i.e.  larger  and  larger  b) 
the  threshold  power  can  be  further  reduced.  However,  one  can  not  reduce 
the  wavelength  to  very  small  values  without  considering  the  collisional 


effects. 
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III.  DECAY  INSTABILITY  IN  PLT 


The  decay  modes  for  a  pump  around  <v  a>c  possess  long  wavelengths 

<  1  (except  for  the  case  of  decay  into  a  Bernstein  mode  and  a  low 

A 

frequency  quasimode  which  has  been  studied  by  Ott  et  al.  in  the  dipole 
approximation)  and  hence,  fluid  theory  can  be  employed  to  obtain  the 
nonlinear  response  of  electrons.  Again,  we  consider  the  propagation 
of  the  extraordinary  pump  wave  [cf.  Eqs.  (1)  and  (2)]  along  x  asix  in 
the  tokamak,  with  a  static  magnetic  field  in  the  z  direction.  The  linear 
response  of  the  pump  may  be  written  as 
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(17) 


The  pump  decays  into  two  electrostatic  waves  <f> (co , k)  and  (j>^(w^,k^) .  The 
response  of  electrons  to  these  fields  is  governed  by  the  equations  of 
motion  and  continuity. 
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Expanding 


n  =  nQ  +  no(u>o,ko)  +  n(w,k)  + 
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the  linear  response  at  w^.k^  may  be  written  as 
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To  obtain  the  nonlinear  response  we  consider  two  cases. 


A)  u)  <  k  v 
'  a,  z  e 


In  this  limit,  the  low  frequency  nonlinearity  arises  through  the 
parallel  ponderomotive  force  F  =  -iek^^. 
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and  hence  the  low  frequency  density  perturbation  comes  out  to  be 
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where  Xg  =  -jz  t1  +  \Tv~  Z*k“v~^  '  EmPloying  Eqs. 
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nonlinear  contribution  to  the  high  frequency  density 


(17)  and  (22),  the 
perturbation  turns  out 


to  be 
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where  *  denotes  the  complex  conjugate.  Using  Eqs.  (22)  and  (23)  in  the 
Poisson's  equation,  we  obtain  Eq.  (13)  as  the  nonlinear  dispersion  relation 
with  the  following  coupling  coefficient. 


V  =  -Xe(l  +  Xt) 
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X^  is  the  ion  susceptibility  at  (w,k). 


(24) 


B)  on  >  k  v 
z  e 


In  this  limit,  u  and  may  be  comparable  and  the  nonlinearity  in 
the  response  at  both  the  frequencies  arises  through  the  equation  of  con¬ 
tinuity  as  well  as  the  ponderomotive  force.  On  writing  equations  similar 
to  Eqs.  (20)  for  the  linear  response  at  (io,it)  evaluating  the  nonlinear  density 
perturbations  and  employing  the  Poisson's  equation,  we  obtain 
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2  111  o  o  o  o 

The  dispersion  relation  Is  again  given  by  Eq.  (13).  Now  we  discuss  some 
specific  channels  of  decay. 

i)  Two  Lower  Hybrid  Waves. 

The  linear  dispersion  relation  for  lower  hybrid  waves,  viz.. 


km.  1/2 


(0  =  7  to..  _  Io>  "  ^(o  ,u>  (1  +  —  o — )] 


2  tou  -  2  iu,u  -u)pltocv 


tells  that  to  <  rain((oc,<op) ,  hence  the  channel  of  decay  involving  two  lower 

hybrid  waves  is  possible  for  2(op  >  w0  ~  wc-  The  decay  waves  can  have  either 

k  >>  k  or  k  <  k  .  In  the  former  case  to  =  so  /2, 
o^o  o 

.2  .  4,4 

k  uu*  to  k 

„  a  _2 _  _2  _z 

U  2  4,4 

u)  a)  k 
o  o 
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and  the  convective  amplification  is 


A  = 


x 


2  2 

where  we  have  assumed  u)Q  'v  (which  corresponds  to  k^/k  =  1/3  for  a  lower 

hybrid  wave  of  w  =  toc/2.  For  PLT  parameters,  =  35  kG,  Ln  10  cm, 

kQ  =  18  cm-1,  k  =  20  kQ,  k/kx  »  10,  A  *  8  x  105u2/c2  «  1.6  at  200  KW/cm2. 

2 

Thus  the  threshold  for  this  decay  is  ^150  KW/cm  . 

For  k  <  k^,  the  decay  wave  frequencies  may  be  considerably  different 

from  each  other.  The  high  frequency  lower  hybrid  wave  (w^,k^)  has  to  be 

2  2 

close  to  u  ,  i.e.,  w  >  ui  »  k,  »  k,  . ;  the  latter  inequality  is  satisfied 
o  p  'V,  c  lz  11 

for  k^  =  k^  =  kQ.  The  amplification  comes  out  to  be 


A  ,  A.  2 

a)  *  k  k  kL 

A  «  _£_  a  o  1  —  kL 

2  2  2  A  8  *  2  n 


0)^.  0). 


2 

This  requires  threshold  powers  above  500  KW/cm  for  PLT  parameters. 

The  decay  instability  of  the  ordinary  pump  wave  possess  almost  the 
same  amplification. 

ii)  Lower  Hybrid  and  Ion  Cyclotron  Waves. 

The  ion  cyclotron  wave. 


w  -  wci(l  + 


-bi 


r+Tj/T 


) 

e 


has  i»>«(o,,u)-u>.>kv.  ,  w  <  k  v  ,  k,  p  <1  and  this  decay  channel  is 
cl  cl  z  i  z  e  l  e 

possible  only  for  an  extraordinary  pump  in  between  -x  <  x  <  0  (where 
2  2 

ui  >  (D  )  and  when  <u  >  w  .  The  resonance  condition  gives  k  >>  k  »  k.. . 

CO  p^o  izli 

i.e.,  k^  *  kQ.  The  convective  amplification  turns  out  to  be 


: 
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A  * 


16v 


,  k  L 
2  y  n 


and  the  threshold  power  >200  KW/cm  for  PLT  parameters, 
lli)  Lower  Hybrid  Mode  and  Quasimode. 

When  the  parallel  phase  velocity  of  the  lower  frequency  lower  hybrid 
wave  equals  the  electron  thermal  speed,  oj  =  ^zve»  it  i-8  strongly  Landau 
damped;  the  higher  frequency  lower  hybrid  wave  is  weakly  damped.  Solving 


3e 


1 


Eq.  (13)  with  a)  ■  w  +  iy,  c,  =  iy  - —  the  homogeneous  medium  growth  rate 

L  3u)  ^ 

turns  out  to  he 


I  (p/e) 

IP _ 

Yo  *  3e1/3o>1 

,  2 
=  3  x  10  (*i 

o  l 
c 

/ 

which  is  'vlO  at  200  KW/cm  and  PLT  parameters.  This  growth  rate  is  extremely 
low  to  substantially  emplify  the  daughter  wave  in  one  pass  through  the  region 
of  parametric  resonance. 

iv)  Filamentation  Instability. 

The  analysis  developed  here  is  applicable  only  to  three  wave  decay 
processes.  However,  one  could  initially  include  the  upper  sideband 
(u>  +  a>0.  t  +t  )  in  the  analysis  to  study  four  wave  processes.  The  fila¬ 
mentation  instability  of  the  ordinary  mode  in  tokamak,  with  the  perturbation 
along  z  axis  (1c  |j  z),  is  not  to  be  affected  by  the  magnetic  field  and  the 
amplification  factor  turns  out  to  be 


A  ^  —  kL 


U)  v 
o  e 


which  again  requires  threshold  power  <v200  KW/cm  for  PLT  parameters.  The 
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fllamentatlon  and  modulation  instabilities  of  the  extraordinary  mode  require 
much  higher  power  densities. 


_ 1 
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IV.  DISCUSSION 


2 

The  presently  employed  microwave  power  densities  (>  200  Watt/cm  ) 
in  EBT  seem  to  be  sufficient  to  parametrically  excite  two  electron  Bernstein 
modes  at  the  second  harmonic  cyclotron  layer,  as  long  as  <  200  eV.  The 
Bernstein  modes  are  strongly  absorbed  through  cyclotron  absorption  and  may 
account  for  the  formation  of  the  hot  electron  annulus.  The  coupling  co¬ 
efficient  for  this  decay  vanishes  in  the  dipole  approximation. 

The  situation  in  large  tokamaks  where  T  >  1  keV  and  B  >  30  kG  is  quite 
unfavorable  for  parametric  instabilities.  The  most  prominant  channel  of 
decay  seems  to  be  the  one  involving  two  lower  hybrid  waves  of  short  wave¬ 
lengths  (k»kQ,  to  -  wo/2),  however,  the  convective  threshold  for  this 

2 

channel  is  ^150  KW/cm  which  Is  orders  of  magnitude  higher  than  the  presently 


employed  power  densities. 
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